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Using the data accumulated in 2002-2004 with the D0 detector in proton-antiproton collisions at 
the Fermilab Tevatron collider with centre-of-mass energy 1.96 TeV, the branching fractions of the 
decays B — > D\ (2420)/i + f M A r and B — ♦ _D2°(2460)^i + ^ M X and their ratio have been measured: 
B(b -> B) ■ B(B -> Dln+v^X) ■ B{&1 -> D*--k + ) = (0.087 ± 0.007 (stat) ± 0.014 (syst))%; 
B(b -> B)j B(B -> Dffi+v^X) ■ B{D? -> D*-tt+) = (0.035 ± 0.007 (stat) ± 0.008 (syst))%; and 
(B(B -> Dffi+^X) ■ B(Df D*-tt+))/(B(B -> D^+u^X) ■ B{D\ -> D*~ tt + )) = 
0.39 ± 0.09 (stat) ± 0.12 (syst), where the charge conjugated states are always implied. 

PACS numbers: 13.25.Hw,14.40.Lb 



This Letter describes our investigation of the prop- excited states of the D meson that have small decay 
erties of semileptonic decays of B mesons to orbitally widths. In the simplest case, these states consist of 
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a charm quark and a light quark in a state with or- 
bital angular momentum equal to one. In the limit of 
a large charm quark mass m c 3> ^qcd, one doublet of 
states with j = 3/2 (Di, D%) and another doublet with 
j = 1/2 (Dq, D-y) are predicted to exist, where the an- 
gular momentum j is the sum of the light quark spin and 
orbital angular momentum. Conservation of parity and 
angular momentum restricts the final states that are al- 
lowed in the decays of these particles collectively known 
as D** mesons. The states that decay through a D-wave, 
D\ and D^, are expected to have small decay widths, 
O(10 MeV/c 2 ), while the states that decay through an 
S-wave, Dq and D 1: are expected to be broad, 0(100 
MeV/c 2 ). 

The ratio R of the semileptonic branching fractions of 
the B meson to D\ and D 2 : 



R 



B(B -> D* 2 lv) 
B(B -> Di£v) ' 



(1) 



is one of the least model-dependent predictions of Heavy 
Quark Effective Theory (HQET) 1] for these states. This 
ratio is expected to be equal to 1.6 in the infinite charm 
quark mass limit Q, but it can have a lower value once 
0(l/m c ) corrections are taken into account 0,0. To- 
gether with the measurement of the corresponding ratio 
Rtt for the non-leptonic decays B — > D**n, determination 
of R will provide important tests of HQET and factor- 
ization of the non-leptonic decays Q . 

The narrow D** mesons have been previously studied 
by several experiments, most recently at Belle |5J where 
the ratio R^ was measured. The semileptonic decay frac- 
tions of B mesons to D** mesons were reported previ- 
ously by the ARGUS @, CLEO 0, OPAL @, ALEPH 
0, and DELPHI (as preliminary) [ljj collaborations, 
with only the latter measuring the fraction of B — > D^lv 
and the others setting upper limits for this decay mode. 

The data set used for this analysis corresponds to 
w 460 pb^ 1 of integrated luminosity accumulated by the 
D0 detector between April 2002 and September 2004 
in proton-antiproton collisions at the Fermilab Tevatron 
collider at centre-of-mass energy 1.96 TeV. The D0 de- 
tector has a central tracking system consisting of a sili- 
con microstrip tracker (SMT) and a central fiber tracker 
(CFT) 0. Both are located within a 2 T supercon- 
ducting solenoidal magnet and have designs optimized 
for tracking and vertexing for \rj\ < 3 and 77 1 < 2.5 
[l2j | , respectively. The SMT has a six-barrel longitudinal 
structure, each with a set of four layers arranged axially 
around the beam pipe, and interspersed with 16 radial 
disks. Silicon sensors have typical strip pitch of 50 — 150 
/im. The CFT has eight thin coaxial barrels, each sup- 
porting two doublets of overlapping scintillating fibers of 
0.835 mm diameter. The next layer of detection involves 
a preshower constructed of scintillator strips and a liquid- 
argon/uranium calorimeter. An outer muon system, cov- 
ering j 77! < 2, consists of a layer of tracking detectors and 



scintillation trigger counters in front of 1.8 T iron toroids, 
followed by two similar layers after the toroids A 
suite of single-muon online triggers was used to record 
the data set while offline only information from the muon 
and tracking systems was used in this analysis. 

Production of narrow D** mesons in B — > 
D*~ tt + n+v^X decay manifests itself as resonance peaks 
in the D*~tt + ll l| invariant mass spectrum. To perform 
the measurement, the semileptonic branching fractions 
of B mesons to the D** mesons were normalized to the 
B — > D*~ [jl+v^X process. 

Initially, a sample of fi^D candidates was selected 
by requiring a muon with transverse momentum pt^ > 
2 GeV/c and 1 77^ | < 2. D° mesons were reconstructed 
through their decays into K + ir~ . Two tracks with pt > 
0.7 GeV/c and \rj\ < 2 were required to belong to the 
same jet and to form a common D° vertex following the 
procedure described in detail in Ref . 0] . To increase the 
signal yield, the event selections of Ref. 0] were relaxed 
by removing the explicit requirement that the px of the 
D° exceeds 5 GeV/c. In total 216870 ±1280 (stat) n+D° 
candidates were found. 

D*~ candidates were selected through their decays into 
D°ir~ by requiring an additional track with pt > 0.18 
GeV/c and the charge opposite to that of the muon. The 
mass difference AM = M(Kirir) — M(Kir) for all such 
tracks with assigned pion mass is shown in Fig. ^ for 
events with 1.75 < M(Kir) < 1.95 GeV/c 2 . The signal 
was described as the sum of two Gaussian functions and 
the background as the sum of exponential and first-order 
polynomial functions. The total number of D*~ candi- 
dates in the peak is 55450 ± 280 (stat) and is defined as 
the number of signal events in the mass difference window 
between 0.142 and 0.149 GeV/c 2 . 

To select a sample of B — > D*~ fi+v^X decays used 
later both for the signal search and for the normalization, 
B candidates were defined using the fi + and D*~ par- 
ticles. All tracks used for the reconstruction of the B 
candidate had to have at least two SMT and six CFT 
hits. The decay length of the B meson, defined in the 
axial plane [l^ | as the distance between the primary ver- 
tex [l5| and the B meson vertex, was restricted to be less 
than 1 cm, the uncertainty on the B- vertex axial position 
had to be less than 0.5 mm, and the x 2 of the £?-vertex 
fit had to be less than 25 for three degrees of freedom. 
The significance of the decay length in the axial plane 
and the proper decay length of the B meson 0] were 
required to exceed 3.0 and 0.25 mm respectively. The 
significance is defined as the ratio of the decay length to 
its uncertainty. The last selection reduces the cc contam- 
ination in the n + D° sample 15]. After these selections, 
the total number of D*~ candidates in the invariant mass 
difference peak is No- = 31160 ± 230 (stat). 

D** decays can be selected by combining the D* can- 
didates with an additional track with assigned pion mass. 
The track was required to have a charge opposite to that 
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FIG. 1: Distribution of the mass difference M(7Y7r7r)-M(/Y7r) 
for events with 1.75 < M(Ktv) < 1.95 GeV/c 2 . The fit func- 
tion describes the signal as the sum of two Gaussian functions 
and the background as the sum of exponential and first-order 
polynomial functions. The signal contribution is also shown 
separately. The hatched histogram corresponds to the same- 
sign combination of the muon and pion charges. 



of the D*, p T > 0.3 GeV/c, and at least two SMT and 
six CFT hits. Pions from the D** decay can also be 
selected by their topology since the corresponding track 
originates from the B-vertex rather than from the pri- 
mary vertex. The impact parameters (IP) in the axial 
plane with respect to the primary and with respect to 
the B-vertex were determined for each track. In order to 
select tracks belonging to the i3-vertex, the ratio of IP 
significances of the track for the primary and B vertices 
was required to be greater than four and the IP signif- 
icance with respect to the primary vertex was required 
to be greater than one. The IP significance is defined as 
the ratio of the impact parameter to its uncertainty. 

The D** invariant mass distribution after all selections 
is shown in Fig. |2 where the D* mass from the Particle 
Data Group (PDG) [13 has been used as a mass con- 
straint. The observed mass peak can be interpreted as 
two merged narrow D** states, D° and D^ - 

The distribution was fit using a sum of two relativistic 
Breit-Wigner functions Gd x and Go*, corresponding to 
the two narrow D** states with the numbers of events 
Nd 1 and Nr>* respectively, and a second-order polyno- 
mial ft describing the background, see Eq. (2). 



fix) 
Gi 

T(x) 

F^(k,k ) 



N Dl ■ G Dl + N D . ■ G DS + f b (x); (2) 

Mi r(x) 



(x 2 - M[) 2 + MfT(x) 2 

2L + 1 



Resi 



r, 



Mi 



x \ko 
9 + 3(k Q z) 2 + (h Q z) 4 
9 + 3(kz) 2 + (kz) 4 



P (i) (*,*o)i 



FIG. 2: The invariant mass M(D*n). Event selection is de- 
scribed in the text. The points correspond to the D*7r combi- 
nations with opposite charges, and the hatched histogram cor- 
responds to the same-charge combinations. The distribution 
was fit with a sum of two relativistic Breit-Wigner functions 
corresponding to two narrow D** states and a second-order 
polynomial describing the background. The contributions of 
D° and Df to the fit are shown separately. 



In the formulas above, x is the D*tt invariant mass 
and Mi and Ti are the mass and width of the corre- 
sponding resonance. The variables k and ko are the pion 
three-momenta in the D** rest frame when the D** has 
a four- momentum-squared equal to x 2 and M 2 , respec- 
tively. F( 2 '(k,ko) is the Blatt-Weisskopf form factor for 
D-wave (L = 2) decays of D** mesons [l8|, and z = 1.6 
(GeV/c) -1 is a hadron scale corresponding to the case 
of the charm quark. Res^ is the mass resolution func- 
tion described by two Gaussian functions with the pa- 
rameterization determined from Monte Carlo (MC) sim- 
ulations. The second Gaussian describing the resolution 
corresponded to 28% of events and was wider by a factor 
of 2.2 than the first one. The standard D0 simulation 
chain included the evtgen Jljj] generator interfaced to 
PYTHIA [2fj| and followed by full GEANT [2l| modeling 
of the detector response and event reconstruction. The 
MC resolution was scaled up by 20% to account for the 
the difference between the data and the MC, where the 
scaling factor was estimated by comparing the D* mass 
resolution in the data and MC. The mass resolution used 
for the fit (sigma of the first Gaussian in the resolution 
function) was 8.2 MeV/c 2 for D\ and 9.4 MeV/c 2 for D*° 
after the scaling. 

The parameters of the background function were de- 
termined by fitting the distribution of the same-charge 
combinations, fixing the function shape parameters, and 
allowing the overall normalization of the background to 
float in the mass fit. The mass difference between 
and D^ and the widths of and D^ were fixed in the 
fit using their corresponding values from the PDG |l7| : 
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TABLE I: Relative systematic uncertainties on the semilep- 
tonic branching fraction to both narrow states, Bd**', the 
semileptonic branching fractions to D° and to D2 , Bd x and 
B D * ; and their ratio B D * /Bd x ■ 



Source 


Bo** 








B(b -> D—l+vX) 


7% 


7% 


7% 




Mass resolution 


2% 


4% 


5% 


8% 




3% 


11% 


16% 


24% 




2% 


2% 


11% 


13% 


AM 


1% 


3% 


6% 


9% 


MC statistics 


2% 


2% 


2% 


3% 


e_D*» modeling 


3% 


3% 


3% 


2% 


Wide resonance 


5% 


5% 


5% 


2% 


Interference effects 


2% 


2% 


2% 


2% 


D* fit and £>** bkg fit 


4% 


4% 


4% 


1% 


Total uncertainty 


12% 


16% 


24% 


30% 



36.7 ± 2.7 MeV/c 2 , 18.9±|J MeV/c 2 , 23 ± 5 MeV/c 2 . 
The numbers of events in the two narrow states de- 
rived from the fit, N Dl = 467 ± 39 (stat) and N D * = 
176 ± 37 (stat), were used to determine the total num- 
ber of events Nd** = 643 ± 38 (stat), and the ratio 
No* /Nd 1 = 0.378 ±0.086 (stat), where the uncertainties 
take into account correlation between the variables. The 
X 2 of the fit at the minimum is 46.9 for 46 degrees of 
freedom. 

The branching fraction for the decays B — > D** [i + v ^X 
can be determined by normalization to the known value 
of the branching fraction B(b -> D*~£+vX) = (2.75 ± 
0.19)% 01 ■ The following two formulas were used for 
the calculations: 

B(b B) ■ B(B -> D**n+v^X) ■ B(D** -> D*-tt + ) = 

B(b^D*-£+vX)-^-—; 

N D * e_D" 

B(B -> Dfn+^X) ■ B(D* 2 ° -> D*-tt+) _ N^_ 
B(B -> • 6(5? -> #*-tt+) ~ A^ ' erjj " 

JVd»« and Nd* are the numbers of D** and D* candi- 
dates as defined above. The D** notation stands for D\ 
or 7?2° or both of them, (D^Df). e D ** is the efficiency 



to reconstruct the charged pion from the D** decay de- 
termined from the MC and equal to (47.2 ± 1.0 (stat))% 
for the D\ meson and (45.4 ± 1.2 (stat))% for the Df 
meson. Contributions from B s mesons, Af, baryons, 
B — > D)ID X decays and the cc process to the sample 
were found to be small and have been neglected 0] . 

The relative systematic uncertainties on the branch- 
ing fractions and of their ratio are summarized in Ta- 
ble [I] The contribution due to uncertainty in Bib — > 
D*~ t + vX) was determined from the uncertainty on this 
branching fraction. The systematic uncertainty caused 
by the MC mass resolution was estimated by varying 
the resolution by ±20%. Contributions due to limited 
knowledge of the D** masses and widths were computed 
by refitting the mass distribution after varying these pa- 
rameters within their uncertainties. The systematic un- 
certainty due to efficiency modeling accounts for the vari- 
ation caused by a possible mismatch between the pr spec- 
tra of reconstructed particles in the data and MC. 

There are predictions and possibly observations |5j of 
a wide resonance D® with a mass of 2430 MeV/c 2 and a 
width of 380 MeV/c 2 predominantly decaying to D*~ir + . 
This resonance is not apparent in our data, and it was 
not used in the fits. The systematic uncertainty caused 
by a possible contribution of this resonance was evalu- 
ated allowing for another Breit-Wigner function in the 
fit, with the mass and width fixed to the wide resonance 
parameters. 

Any interference effects between the and Df must 
average to zero after integration over all angles under the 
assumption of equal acceptances. The validity of this as- 
sumption has been checked and the corresponding uncer- 
tainty assigned. The systematic uncertainty due to the 
fitting procedure was estimated by varying the functions 
describing the backgrounds for the D* and D** mass dis- 
tributions and also the function describing the D* mass 
peak. The total systematic uncertainty was found by 
summing all the above sources in quadrature. 

Using the numbers defined above, the semileptonic 
branching fractions of B mesons to D** mesons and their 
ratio are: 



B{b -» B) ■ B(B -> (Dj, Df)^i+v^X) ■ B{{D\, Df) -» D*~ir + ) = (0.122 ± 0.007 (stat) ± 0.015 (syst))%; 

B(b -> B) ■ B(B -> Dlfi+v^X) ■ B{D° -> D*-tt+) = (0.087 ± 0.007 (stat) ± 0.014 (syst))%; 

B(b -> B) ■ B(B -> Dfn+v^X) ■ B(D? -> D*~-k+) = (0.035 ± 0.007 (stat) ± 0.008 (syst))%; 
B(B -> Dffx+Vf.X) ■ B{Df -> D*-ir+) 



B(B -» D° lt i+^X) ■ B{Dl -> D*-^ 



0.39 ± 0.09 (stat) ± 0.12 (syst). 



Upon using the input B(b —>£?) = (39.7 ± 1.0)% [J?] , assuming isospin conservation and that the D± me- 
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son decays only into D*ir |22|, the branching fraction 
B(B -> Dll+vX) = (0.33 ± 0.06)% is determined. It 
is different from the PDG value, (0.74 ± 0.16)% [HI, by 
2.5 standard deviations. Similarly the branching frac- 
tion B(B -► Dfl+vX) = (0.44 ± 0.16)% is determined 
assuming that D 2 decays into D*tt in (30 ± 6)% of the 



cases 1YJ . This result is in agreement with the 95% CL 
upper limit 0.65% provided by the PDG. 

Using the measured ratio of the branching fractions 
and the same assumptions for the absolute fractions for 
Di and D* 2 as above, the ratio R — 1.31 ± 0.29 (stat) ± 
0.47 (syst) was computed. 

In summary, using 460 pb _1 of integrated luminosity 
accumulated with the D0 detector, the semileptonic de- 
cays B — > Dlii+v^X and B — > D^ 1 ^l + v^X have been ob- 
served and the branching fractions measured using statis- 
tics more than an order of magnitute better than previ- 
ous measurements 0, S El • This result represents 
a significant improvement in the knowledge of B branch- 
ing fractions to orbitally excited D mesons, and the first 
direct measurement of their ratio. 
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